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We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at z > 0.3, 
which cover a temperature range of 3 > kT > 15 keV. Our analysis is aimed at measuring 
the iron abundance in the ICM out to the highest redshift probed to date. We find that the 
emission- weighted iron abundance measured within (0.15 — 0.3) R V i r in clusters below 5 keV 
is, on average, a factor of ~ 2 higher than in hotter clusters, following Z(T) ~ 0.88 J 1-0 - 47 Zq, 
which confirms the trend seen in local samples. We made use of combined spectral analysis 
performed over five redshift bins at 0.3 > z > 1.3 to estimate the average emission weighted 
iron abundance. We find a constant average iron abundance Zfb — 0.25 Zq as a function 
of redshift, but only for clusters at z > 0.5. The emission- weighted iron abundance is 
significantly higher (Zfc — 0.4 Zq) in the redshift range z ~ 0.3 — 0.5, approaching the value 
measured locally in the inner 0.15 R v i r radii for a mix of cool-core and non cool-core clusters 
in the redshift range 0.1 < z < 0.3. The decrease in Zp e with z can be parametrized by 
a power law of the form ~ (1 + z) 125 . The observed evolution implies that the average 
iron content of the ICM at the present epoch is a factor of ~ 2 larger than at z ~ 1.2. We 
confirm that the ICM is already significantly enriched [Zfh ~ 0.25 Zq) at a look-back time 
of 9 Gyr. Our data provide significant constraints on the time scales and physical processes 
that drive the chemical enrichment of the ICM. 
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§1. Properties of the sample and spectral analysis 



^ ■ The selected sample consists of all the public Chandra archived observations of 

clusters with z > 0.4 as of June 2004, including 9 clusters with 0.3 < z < 0.4. We 
used the XMM.- Newton data to boost the S/N only for the most distant clusters in 
our current sample, namely the clusters at z > 1. 

The spectrum of each cluster is extracted from a circular region whose radius 
maximizes the S/N. As shown in Fig. [IJ in most cases the extraction radius R ex t is 
between 0.15 and 0.3 R v i r . The spectra were fitted with a single-temperature model 
in which the ratio between the elements was fixed to the solar valueP 

We show in Fig. Q] the distribution of temperatures in our sample as a function 
of redshifts (error bars are at the la c.L). The Spearman test shows no correlation 
between temperature and redshift (r s = —0.095 for 54 d.o.f., probability of null 
correlation p = 0.48). Fig.[T]shows that the range of temperatures in each redshift bin 
is about 6 — 7 keV. Therefore, we are sampling a population of medium-hot clusters 
uniformly with z, with the hottest clusters preferentially in the range 0.4 < z < 0.6. 

Our analysis suggests higher iron abundances at lower temperatures in all the 
redshift bins. This trend is somewhat blurred by the large scatter. We find a more 
than 2a negative correlation for the whole sample, with r s = —0.31 for 54 d.o.f. 
(p = 0.018). The correlation is more evident when we compute the weighted average 
of Zp e in 6 temperature intervals, as shown by the shaded areas in Fig. [2J 
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Fig. 1. Left: extraction radius R ex t versus R v i r . Lower and upper lines show R ex t = 0.15 R V i r 
and Rext = 0.3 R V i r , respectively. Right: temperature vs redshift. Shaded areas show the rms 
dispersion around the weighted mean in different redshift bins. 



§2. The evolution of the iron abundance with redshift 

The single-cluster best-fit values of Zp e decrease with redshift. We find a ~ 3a 
negative correlation between Zp e and z, with r s = —0.40 for 54 d.o.f. (p = 0.0023). 
The decrease in Zp e with z becomes more evident by computing the average iron 
abundance as determined by a combined spectral fit in a given redshift bin. We 
performed a simultaneous spectral fit leaving temperatures and normalizations free 
to vary, but using a single metallicity for the clusters in a narrow z range. 

The Zp e measured from the combined fits in 6 redshift bins is shown in Fig. [2 We 
also computed the weighted average from the single cluster fits in the same redshift 
bins. The best-fit values resulting from the combined fits are always consistent with 
the weighted means within la (see Fig. [2|). This allows us to measure the evolution 
of the average Zp e as a function of redshift, which can be modelled with a power 
law of the form ~ (1 + z)~ L25 . 

Since the extrapolation of the average Zp e at low-z points towards Zp e (0) ~ 
0.5 Zq, we need to explain the apparent discrepancy with the oft-quoted canonical 
value (Zp e ) ~ 0.3 Zq. The discrepancy is due to the fact that our average values 
are computed within r ~ 0.15 R V i r , where the iron abundance is boosted by the 
presence of metallicity peaks often associated to cool cores. The regions chosen for 
our spectral analysis, are larger than the typical size of the cool cores, but smaller 
than the typical regions adopted in studies of local samples. In order to take into 
account aperture effects, we selected a small subsample of 9 clusters at redshift 
0.1 < z < 0.3, including 7 cool-core and 2 non cool-core clusters, a mix that is 
representative of the low-z population. Here we analyze the X-ray emission within 
r = 0.15i?„i r in order to probe the same regions probed at high redshift. We used 
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Fig. 2. Left: scatter plot of best-fit Zf b values versus kT. The dashed line represents the best-fit 
Z-T relation (Z/Z Q ~ 0.88 T -0 ' 47 ). Shaded areas show the weighted mean (blue) and average 
ZFe with rrns dispersion (cyan) in 6 temperature bins. Right: mean Zfb from combined fits (red 
circles) and weighted average of single-source measurements (black squares) within 6 redshift 
bins. The triangles at z ~ 0.2 are based on the low-z sample described in Sect. 2. Error bars 
refer to the la c.L. Shaded areas show the rms dispersion. The dashed line indicates the best 
fit over the 6 redshift bins for a simple power law of the form (Zj? e ) = Zp e (0) (1 + z)~ 1,25 . 



this small control sample to add a low-z point in our Fig. [21 which extends the Zp e 
evolutionary trend. 

§3. Conclusions 

We have presented the spectral analysis of 56 clusters of galaxies at intermediate- 
to-high redshifts observed by Chandra and XMM- Newton^ This work improves our 
first analysis aimed at tracing the evolution of the iron content of the ICM out to 
z > l, 4 ' by substantially extending the sample. The main results of our work can 
be summarized as follows: 

• We determine the average ICM iron abundance with a ~ 20% uncertainty at 
z > 1 {Zp e = 0.27 ± 0.05 Zq), thus confirming the presence of a significant 
amount of iron in high-z clusters. Zp e is constant above z — 0.5, the largest 
variations being measured at lower redshifts. 

• We find a significantly higher average iron abundance in clusters with kT < 
5 keV, in agreement with trends measured in local samples. For kT > 3 keV, 
Zp e scales with temperature as Zp e (T) ~ 0.88 T~ 0A7 . 

• We find significant evidence of a decrease in Zp e as a function of redshift, 
which can be parametrized by a power law {Zp e ) — -^Fe(O) (1 + z)~ az , with 
Zp e (0) ~ 0.54 ± 0.04 and a z ~ 1.25 ± 0.15. This implies an evolution of more 
than a factor of 2 from z = 0.4 to z = 1.3. 
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We carefully checked that the extrapolation towards z ~ 0.2 of the measured 
trend, pointing to Zp e ~ 0.5 Z®, is consistent with the values measured within a 
radius r = 0.15 Rmr i n local samples including a mix of cool-core and non cool- 
core clusters. We also investigated whether the observed evolution is driven by a 
negative evolution in the occurrence of cool-core clusters with strong metallicity 
gradients towards the center, but we do not find any clear evidence of this effect. We 
note, however, that a proper investigation of the thermal and chemical properties of 
the central regions of high-z clusters is necessary to confirm whether the observed 
evolution by a factor of ~ 2 between z = 0.4 and z = 1.3 is due entirely to physical 
processes associated with the production and release of iron into the ICM, or partially 
associated with a redistribution of metals connected to the evolution of cool cores. 

Precise measurements of the metal content of clusters over large look-back times 
provide a useful fossil record for the past star formation history of cluster baryons. 
A significant iron abundance in the ICM up to z ~ 1.2 is consistent with a peak in 
star formation for proto-cluster regions occurring at redshift z — 4 — 5. On the other 
hand, a positive evolution of Zp e with cosmic time in the last 5 Gyrs is expected on 
the basis of the observed cosmic star formation rate for a set of chemical enrichment 
models. Present constraints on the rates of SNae type la and core-collapse provide a 
total metal produciton in a typical X-ray galaxy cluster that well reproduce (i) the 
overall iron mass, (ii) the observed local abundance ratios, and (iii) the measured 
negative evolution in Zp e up to z ~ 1.2P 
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